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Numerical integration of a stochastic Landau-Lifshitz-Gilbert equation is used to study dynamic 
processes in single-domain nanoscale magnets at nonzero temperatures. Special attention is given 
to including thermal fluctuations as a Langevin term, and the Fast Multipole Method is used to 
calculate dipole-dipole interactions. It is feasible to simulate these dynamics on the nanosecond time 
scale for spatial discretizations that involve on the order of 10* nodes using a desktop workstation. 
The nanoscale magnets considered here are single pillars with large aspect ratio. Hysteresis-loop 
simulations are employed to study the stable and metastable configurations of the magnetization. 
Each pillar has magnetic end caps. In a time-dependent field the magnetization of the pillars is 
observed to reverse via nucleation, propagation, and coalescence of the end caps. In particular, 
the end caps propagate into the magnet and meet near the middle. A relatively long-lived defect is 
formed when end caps with opposite vorticity meet. Fluctuations are more important in the reversal 
of the magnetization for fields weaker than the zero-temperature coercive field, where the reversal is 
thermally activated. In this case, the process must be described by its statistical properties, such as 
the distribution of switching times, averaged over a large number of independent thermal histories. 



The effect of temperature on the switching behav- 
ior of nanoscale magnets can be quite strong when ex- 
ternal fields are applied that are just below the zero- 
temperature coercive threshold. Under these conditions, 
thermal fluctuations can provide enough energy to take 
the magnetization of the system over the barrier that pre- 
vents it from aligning with the external field 0. These 
issues are important for understanding data integrity and 
high-speed switching in single-domain magnetic applica- 
tions. 

Nanoscale magnets are modeled using the traditional 
Landau-Lifshitz-Gilbert equation 

whereby the microscopic dipoles Mi precess under their 
individual, locally observed applied fields H^. The uni- 
versal gyromagnetic factor is 7o=1.76 x 10'' Hz/Oe, 
while the material parameters were selected to match 
those of bulk iron with a saturation magnetization 
Ms=1700 emu/cm'^, exchange length 1^=2.6 nm, and 
damping parameter a=0.1. For the results considered 
here, these fields are composed of contributions from a 
uniform field external to the system, the exchange inter- 
action with neighboring dipoles, and dipole-dipole inter- 
actions with all of the other dipoles in the system. Cal- 
culation of the latter is dramatically accelerated by using 
a fast multipole algorithm ||^ . The numerical models are 
based on real nanomagnets that have been fabricated re- 
cently using scanning microscopy techniques Q . 

Thermal effects are incorporated by adding a ran- 
dom contribution to the local field of each spin, as first 
proposed by W. F. Brown ^ nearly forty years ago. 



He calculated the fluctuation-dissipation theorem for a 
Stratonovich-type 1^ noise by considering the Fokker- 
Planck formulation of the problem. Numerically we 
implemented this noise using an Ito-type random 
field. The difference between the Stratonovich and Ito 
paradigms does not affect the results since we normalize 
the dipoles to a flxed length after each integration step. 

The individual magnets, rectangular with dimensions 
9 nm X 9 nm x 150 nm, were implemented using a finite- 
difference approach. The numerical results have been 
found to be independent of discretization for the spa- 
tial discretization of Aa;=1.5 nm and integration step of 
Ai=50 fs used here. 

The large shape anisotropy of these magnets causes 
spontaneous alignment of the magnetization, except at 
the top and bottom where pole avoidance leads to the for- 
mation of end caps. As can be seen in Fig. 1(a), when the 
external field is applied in the opposite direction during 
a hysteresis measurement, the z-component of the mag- 
netization is reduced via growth of these end caps. Here 
the uniform applied field points down, light shades in- 
dicate upward-pointing magnetization, and dark shades 
indicate downward-pointing magnetization. At zero tem- 
perature the end caps grow symmetrically; when they 
meet at the middle of the pillar a relatively long-lived de- 
fect is formed due to the opposite helicities of the two end 
caps. These results are fully consistent with the T = 
simulations of Ref. . There the field was swept quasi- 
statically as opposed to the truly dynamic sweeps con- 
sidered here. 

The corresponding hysteresis loop, with a period of 
1 ns, is shown in Fig. 2. The defect is indicated by slow 
decay of the magnetization around the time when the ex- 
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ternal field reverses. The effect of this defect is also ap- 
parent for the longer-period hysteresis curves, for which 
the defect disappears before the field reverses. The sim- 
ulated zero-temperature hysteresis loops are very repro- 
ducible, and no differences are seen for subsequent peri- 
ods. The inset shows example variations that occur from 
thermal effects for the 1 ns hysteresis loop. The largest 
differences are seen in conjunction with the defect that 
forms from the two end caps. 

The simplest magnetization-reversal situation to study 
is one in which the external field suddenly changes its 
orientation, and then remains constant. In these simula- 
tions, the field is initially zero and then is brought nega- 
tive in 0.25 ns with its amplitude described by 1/4 of a 
sine wave. In what follows we set t=0 at the time when 
the field first reaches its maximum negative value. When 
the final field is less than the coercive field, the magne- 
tization remains oriented upward until thermal fiuctua- 
tions take it over the associated free-energy barrier. In 
long pillars the end caps do not interact strongly, and the 
free energy of each can be considered separately. The 
free energy as a function of end cap volume has essen- 
tially three cxtrcma: one local minimum corresponding 
to a small end cap, one local maximum corresponding to 
an unstable volume where the tendency for shrinkage is 
equal to the tendency for growth, and the global mini- 
mum corresponding to a switched magnet (spanning end 
cap). Snapshots from the magnetization reversal of one 



Taking p and v as parameters, a nonlinear fit of this two- 
exponential decay theory to the first two moments of the 
data is shown as the dashed line in Fig 3. For compari- 
son, a similar fitting has also been performed for an error 
function form with two parameters (corresponding to a 
Gaussian histogram of switching times), which is shown 
as the dotted curve in Fig. 3. Prom the 85 switches pre- 
sented here there is no clear advantage to either fitting 
function. Different combinations of field and tempera- 
ture should probe regions where the fitting functions are 
not so similar. 

In summary, Landau-Lifshitz-Gilbert dynamics have 
been simulated for three-dimensional models of single- 
domain nanoscale magnets of large aspect ratio. 
Hysteresis-loop and field-reversal simulations show that 
magnetization reversal occurs through the nucleation, 
growth, and coalescence of the end caps. For field- 
reversal simulations at nonzero T which require thermal 
fluctuations to complete the reversal, a simple theory 
that considers the nucleation rate and growth velocity 
of the end caps adequately describes the statistical dis- 
tribution of the switching times. 
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magnet in this situation at _ff=— 1850 Oe and 20 K are 
shown in Fig. 1(b). Here the lower end cap undergoes a 
large fluctuation to take it past the critical volume, after 
which it grows at an almost constant rate to fill the entire 
magnet. (Again a long-lived defect forms when the two 
end caps come into contact.) 

The thermally-activated nature of the end cap growth 
leads to a distribution in the switching times, defined as 
the time when Mz=0. The probability of not switch- 
ing, P„ot(t), for 85 switches at H = -1800 Oe and 20 K 
is shown as the heavy, stepped curve in Fig. 3. Under 
these conditions, the majority of simulated switches oc- 
cur between 0.5 and 1.2 ns after the field reversal finishes. 
There are no switches before 0.4 ns because it takes this 
amount of time for a single supercritical end cap to grow 
to fill half of the magnet. The traces of the average pillar 
magnetization density as a function of time are shown for 
five different switches in the inset of Fig. 3. There are es- 
sentially two slopes observed during the actual switching 
process, corresponding to cases where one or both end 
caps are growing. 

The observation that the end caps decay essentially 
independently and exponentially, with rate p, and that 
freely growing end caps change the global magnetization 
(normalized to lie between ±1) at a constant rate v, can 
be used to construct a simple model to describe the dis- 
tribution of switching times observed in individual exper- 
iments. The resulting probability of not switching is 

(2) 
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FIG. 1. Visualization of the z-component of the magnetization for 9 nm x 9 nm x 150 nm iron nanomagnets (shown in 
3/4 cut-away view). Light shades indicate upward-pointing magnetization, while dark shades indicate downward-pointing 
magnetization, (a) Switching via symmetric growth of end caps for 1 ns period hysteresis loop at K. (b) Switching at 
20 K and ii"=1850 Oe via thermal flu ctuations of one end cap over the saddle point, with subsequen t growth to switch the 
magnetization. Movies are available at tittp : //www. scri . fsu. edu/~browgnrg/micromag/pillar .html 
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FIG. 2. Simulated hysteresis loops at K for applied fields 
oscillating with periods of 1, 2, and 4 ns. The inset shows 
variations caused by temperature for two loops simulated at 
20-?^. The feature that occurs at nearly saturated magnetiza- 
tions is due to a defect that forms when the end caps come in 
contact. 




tsw (ns) 

FIG. 3. Probability of not switching, Pnot, for nano- 
magnets experiencing external fields just below the coercive 
value. The solid line is simulation data from 85 switches at 
ff=-1800 Oe and 20 K. Two theoretical forms [the dotted 
line is an error function corresponding to a Gaussian his- 
togram of switching times, and the dashed line corresponds 
to the two-exponential decay theory described in the text, 
Eq. (2)] are compared after nonlinear fitting with two param- 
eters. The inset shows 5 examples of the average z-component 
of the magnetization during different switching events. 
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